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the layers, mediated by the metal atom. The bond energies between the various atom pairs are 
analysed in terms of orbital- and energy-resolved contributions. 
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I. INTRODUCTION 

The ceramic compound magnesium diboride, MgB2, 
has the highest superconducting critical temperature T c 
= 39 K ever reported for a binary system. The structure 
is graphite-like, i.e., it consists of honeycomb layers q£ 
B separated by triangular metal planes. Its discoveryEJ 
has triggered a tremendous ^ctiszity to investigate the 
electronical] and phononicH'HflO-properties as well as 
the electron-phonon couplingauEnlj for MgB2 and for re- 
lated binary and ternary borides. The main electronic 
characteristics are: 

• The states at the Fermi level Ep are dominantly de- 
rived from boron p orbitals. The B p x>y cr-bands are 
responsible for the strong covalent bonding in the 
B layers and have mainly 2D character, i.e., only 
very little k z dispersion, whereas the B p z 7r-bands 
have 3D character. 

• Electrons are transferred from the Mg orbitals to 
the B orbitals. As a result, there is a partial ionic 
character, but these electrons are not well localized 
at the B sites, they are rather distributed over the 
whole crystal. 

• The Mg atoms are responsible for a downward shift 
of the B 7r-bands relative to the cr-bands, as com- 
pared to the corresponding C bands in graphite. 
As a result, there is a partial occupation of the 7r- 
bands and a hole doping of the cr-bands. Coupling 
of the cr-holes to the intralayer B bond-stretching 
phonon modes dominates the electron -nhonon cou- 
pling and is responsible for the high T c cm. 

• The compound is hold together by dominant in- 
tralayer B-B covalent bonds, dominant interlayer 
metallic bonds and by a substantial ionic charac- 
ter. There are indications for an additional inter- 
layer B-Mg contribution to the bonding, e.g., the 
slight increase in the k z dispersion of the B cr-bands 
as compared to a hypothetical structure where the 



Mg atoms are removed. When going from MgB2 
to AIB2, the valence charge-diference plots show a 
stronger directional M-B bondO (M=metal atom), 
and our own calculations show a slightly stronger 
k z dispersion of the cr-bands. 

It is the objective of the present paper to analyse and 
compare the covalent bonding properties of MgB2 and 
AIB2 quantitatively within a variant, pf our recently de- 
veloped energy-partitioning schemdi3 : Lj. A first attempt 
in this direction was made by Ravindran ef-aZ.El by means 
of the crystal-orbital-Hamilton populatiorJla (COHP). It 
will be shown below, however, that for the use of non- 
orthogonal basis functions the COHP do not have a well- 
defined physical meaning and cannot be used to compare 
quantitatively the covalent bonding properties between 
various structures. 



II. THE ENERGY-PARTITIONING SCHEME 

Starting point of the discussion is the expression for 
the total energy E to t from the density functional theory 
in local density approximation, 

Etot = ^2 f n£n ~ / ™( r Hff(r)d 3 r 
n J 

+ E n + £xc + J n(r)v ext (r)d 3 r + E h , (1) 

with the occupation numbers f n and the eigenvalues e n 
of the Kohn-Sham single-particle wave functions the 
electron density n(r), the effective potential v e s(r), the 
Hartree energy £"h> the exchange-correlation energy E xc , 
the potential v ex t(r) of the nuclei (or of the ionic cores 
in the case of a pseudopotential calculation) and the in- 
teraction energy En between the nuclei (or between the 
ionic cores). The hope is that the trends for the to- 
tal energies are already well described by the first term 
(band-structure energy i^band) when comparing various 
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systems. This assumption is made implicitly in the com- 
mon practice to discuss the energetics via the electronic 
density of states. An energy-partitioning scheme is a tool 
to analyse -Eband further. 

For systems with covalent bonding the crystal wave- 
functions ^ n are best represented by a minimal set of 
well-localized orbitals (fi a (with character a described 
by the angular and magnetic atom quantum numbers I 
and m) attached to the various atoms i (tight-binding 
representation) , 



■•ph 



(2) 



The band-structure energy is calculated by inserting (g) 
into 

E/„£ n =E/n(*n|H|* n > (3) 



with the Kohn-Sham Hamiltonian H. The remaining 
terms on the right-hand side of eq. (Q) which we abbre- 
viated by D in Refs. [l3|,[l4| are calculated by approximat- 
ing n(r\— by a superposition of atomic charge densities, 
yieldingpJ 



D = R 



pair 



E, 



mb 



E 



free atom 
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(4) 
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This term describes the cost in energy when start- 
ing the condensation process from free atoms and then 
redistributing the electrons among the various orbitals 
from the occupation numbers N[™ e atom to the occupa- 
tion miHiber qi a found in the crystal and characterized 
by gio-ta. The second term is the crystal-field term 



Erf — ^ ^ Qia (^Eli a i a Ef r< 



free atom\ 



(9) 



which describes the change in energy due to a shift of the 
on-site energies when the atoms are condensed to form 
the crystal so that the potential acting on an electron 
at atom i is not just the atomic potential of this atom 
but the environment-dependent crystal potential. The 
polarization energy 

-Epolar = E fn C™ Q (c-Jj) Hiai0 ~ 5\^H iaia (10) 

describes the change in energy due to the hybridization of 
orbitals localized at one atom when the atom is embed- 
ded in the crystal. Finally, the covalent bond energy E cov 
is the change in energy arising from the hybridization of 
orbitals localized at different atoms, 



Here -Epair is a pair potential term, E^ a many-body 
potential term which is small for nearly charge neutral 
systems, E hcc atom denotes the total energy of the free 
atoms before being condensed to the crystal, and the 



N^ e atom are the occupation numbers for the energy lev- 

„l c rrfrcc atom 
elb n %aia 

terms 



We now add to and subtract from Rot the 



free atom^ 
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n ia.J/3 



j/3,ia 



H 



free atom 



la/a i 
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where Ojpi a and Hjpi a denote the elements of the over- 
lap and Hamiltonian matrix. 

Rearranging the terms in an appropriate manner yields 
the cohesive energy 

E c = E tot -E trccatom (6) 

— -Eprom H~ -t^cf ~t~ ^polar H~ ^cov H~ -^pair ^mb 

The first term is the promotion energy 



rp \ ^ / Arfree atom\ rrfree atom 

^prom — / 4 \<lia iV i a ) ^iaia 



(7) 



with Mulliken's gross charge 
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cov ,iaj{3 



(ii) 



with 



E 



cov, iotjfi 



E f n C *« (tfp) \- H iP ia ~ OjPia £j/3ia] (12) 



(13) 



E CO v,iajf3 can be subdivided further into energy-resolved 
contributions, 



E cav , iaj p{E) = E S ( E - £„) /„ c" a (c"^) 

n 



(14) 



E CO v,iajf3(E) is negative (positive) for bonding (antibond- 
ing) states. The respective quantity integrated up to a 
certain energy E will be referred to as 



(E>) 

J — CO 



AE' 



(15) 



The energy-partitioning scheme discussed above has 
the following very attractive property: In a band- 
structure calculation which deals with an infinitely ex- 
tended periodic system the average effective potential 
does not have a physical meaning, and it is therefore 
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set to an arbitrary value which is the same for different 
crystal structures. In order to be physically meaning- 
ful the total energy E to t and the considered terms of an 
energy-partitioning scheme for a band-structure calcula- 
tion therefore must be invariant against a constant shift 
of the effective potential of the band-structure calcula- 
tion. This is fulfilled for E to t (and hence also for E c ) 
which becomes obvious from eq. (Jl|) : Shifting v c g by $o 
yields opposite shifts for the first two terms which there- 
fore compensate each other (the remaining terms can be 
calculated without ambiguity, see for instance, Ref. Rh. 
Furthermore the terms Ep rom , E po i ar and E cov of eq. (g) 
as well as their atom- and orbital-resolved contributions 
(and in addition the respective energy-resolved contribu- 
tions to -Ecov) are all invariant. For instance, if the poten- 
tial is shifted by a constant <&0j then the matrix elements 
Hjflia are transformed into Hjpi a + §oOjpi a and Hi a i a , 
£j/3ia into H iaia + $o, + $o because O ia i a = 1, so 

that $o drops out of the covalent bond energy. Because 
E c is also invariant, this must hold for the sum of the 
terms E pair + E m b + E c f, too. However, we cannot cal- 
culate separately, for instance, the crystal-field term E c f 
in a band-structure calculation, because the matrix ele- 
ment Hiaia is shifted by the shift of the effective potential 
of the crystal whereas Hl a c fa tom is not (because for the 
calculation of the latter quantity the effective potential 
is always normalized to zero for distances far from the 
atom). It is therefore physically meaningful only to dis- 
cuss the terms E prom , E po i ar , E cov and E pair + E mb + E c f. 
The covalent bond energy thereby is the only term which 
involves matrix elements for orbitals on different atoms, 
and therefore it clearly represents the contribution of the 
interatomic bonding to the band-structure energy. In 
the following we will confine ourselves to the discussion 
of this term. 

It should be noted that iuij-Qur former version of the 
energy-partitioning schemetJtj the various terms haiie 
been arranged in a slightly different manner, arrivingE3 
at the equivalent expression 

E c = E prom + E c f + E cov + E pa i r + E m b , (16) 
with 

Ep r0 m — ^ ^ (^Qia ^ia ) Eliaia ; (17) 

ia 

i.e., i?£™ atom of E proni has been replaced by H iaia , 

ttt \ A AT-frcc atom / tt rjfrcc atom\ /io\ 

Ecf - N ia {Hiaia ~ H. lma ) , (18) 

ia 

i.e., the qt a of E c f have been replaced by Nf a cc atom ; and 

E CO v — ^ E cov iaj f3 i (19) 

ia,j/3 

i.e., the on-site contributions iai(3 have not been excluded 
from the covalent bond energy. When comparing the 



definitions E prom , E c f, with the definitions E prom , E c f, 
it becomes obvious that they correspond to a different 
succession of processes in a gedanken experiment for the 
condensation of free atoms to the crystal. In the first case 
we promote the electrons by redistributing them among 
the various orbitals of the free atoms and then we bring 
the free atoms in the crystal positions (without allowing 
for a redistribution of the charge densities) and experi- 
ence a change in energy described by E c f due to a shift of 
the on-site energies in the crystal potential. In the second 
case we freeze the occupation numbers N^ a cc atom when 
bringing the free atoms to the crystal positions and cal- 
culate the crystal-field shift E c f for these circumstances, 
and then we allow for a redistribution of the electrons 
among the on-site energy levels in the crystal potential. 
We think that the first case is closer to the commonly 
used definitions of the promotion and the crystal-field 
energy, and we therefore prefer the new variant of the 
energy partitioning scheme. Finally, we think that it is 
reasonable to exclude the on-site hybridization contri- 
butions E cov ,i Q ;^ from the covalent bond energy because 
they do not describe interatomic interactions. In the new 
variant these terms enter the polarization energy which 
also has a well-defined physical meaning. 

It should be noted that E prom , E cov and E c f + E pa ; r + 
EJmb are also invariant against a constant shift of the ef- 
fective potential of a band-structure calculation, as well 
as the atom-, orbital- and energy-resolved contributions 
of E cov . Again, E c f alone is not invariant and cannot be 
calculated separately in a band-structure calculation. 

In former papers the crystal-orbital-Hamilton 
populationta, 

C0HP ia]P {E) =J2 S ( E ~ £ n)/«C {ch)* H m* , (20) 

n 

has been used to characterize the bonding properties. 
If orthonormal basis functions are used, this quantity is 
identical to E cm ; la ip{E) . However, in the chemical anal- 
ysis very often non-orthogonal basis sets are used. Then 
COHPiajp(E) is not invariant against a constant shift 
of the effective potential and therefore does not have a 
well-defined physical meaning in the context of a band- 
structure calculation. 

For systems with metallic bonding the ^> n are better 
represented by a set of plane waves rather than by atom- 
localized functions. Alternatively, the ^f n again can be 
represented by a set of atom-localized functions also in 
this case, but then often orbitals have to be included 
which are not occupied in the free atom in order to make 
the basis set more complete. Then formally a covalent 
bond energy E cov can be calculated even for a nearly-free- 
electron system. Whereas the term metallic conductivity 
is well defined, it is indeed a problem to discriminate 
between metallic and covalent bonding. As a working 
hypothesis we define a covalent bonding as a bonding 
which is dominated by the hybridization of those orbitals 
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on various atoms which are already occupied in the re- 
spective free atoms. Note that this does not necessar- 
ily mean that the corresponding charge-density difference 
plot exhibits directionality, e.g., it can be imagined that 
the p z -orbitals do not necessarily create a charge density 
with considerable directionality. Because our definition 
of -E C ov is a generalization-of the covalent bond energy in- 
troduced by Sutton et alS3 to the case of non-orthogonal 
basis sets we keep this historical founded nomenclature 
"covalent bond energy" although this quantity may also 
contain metallic bonding contributions in the above de- 
fined sense. 



III. DETAILS OF THE CALCULATIONS 

The calculations were-, performed with the ab-initio 
pseudopotential methocE2rE3. Band structure and charge 
density difference plots wexe practically identical to those 
of previous calculations] U. Fqe the energy-partitioning 
analysis the ^ n were projectecOLJ onto a set of overlap- 
ping atom-localized non-orthogonal orbitals. For these 
orbitals we chose 



Via(r) = fu(r)i l K Zm (r) 



(21) 



(1 



-7«( r 



r ) 2 ) for 



r < rT* 
for r > r|; ut , 



(22) 

where Cn is a normalization constant, <^ s is the ra- 
dial pseudo-atomic wavefunction constructed according 
to VanderbiltEU, Xa denotes a contraction factor and r™* 
represents a cut-off length. The parameters Xu, y# 
r^ ut were selected in such a way that the spillage 
was minimized, where the spillage characterizes the loss 
of the norm of the wavefunctions due to the incomplete- 
ness of the pseudo-atomic-orbital projection. We con- 
fined ourselves to a minimal set of s,p and d orbitals for 
Mg and B. The band-structure calculated with the pro- 
jected wavefunctions was nearly identical to the original 
band structure from the pseudopotential calculation for 
energies below and not too far above the Fermi level. 



IV. RESULTS AND DISCUSSION 

Table || represents the covalent bond energies for vari- 
ous atom pairs in MgB2 and AIB2. It should be recalled 
that the covalent bonding properties of the two materi- 
als may be compared only by the measure E cov and not 
by COHP because the latter quantity is not invariant 
against a constant shift of the potential when nonorthog- 
onal basis functions are used, which is the usual case for 
a chemical analysis in terms of atom-like functions. The 
total covalent bond energy is larger in absolute value for 



AIB2 than for MgB2. As expected, the B-B intralayer 
covalent bond energy is largest. Interestingly enough, it 
is larger in absolute value for the Mg compound although 
the B-B distance is larger. Contrary to the general belief, 
there is also a considerable covalent bond energy between 
the layers (the B-M energy is only a factor of about 2-2.5 
smaller) and even in the M layers (the intralayer M-M 
energy is a factor of about two smaller than the B-M 
energy). The B-M covalent bond energy increases when 
going from M=Mg to M=A1, and this is consistent with 
the smaller c/a ratio of the Al compound. The B-M dis- 
tance is in fact smaller for M=A1 than for M=Mg. The 
interlayer B-B (M-M) bonds are two orders of magnitude 
(one order of magnitude) smaller than the respective in- 
tralayer bonds, i.e., their contributions to the cohesion 
between the layers is much smaller than the B-M con- 
tribution. Our calculations show that the covalent bond 
energies between all further-distant atoms are consider- 
ably smaller than the nearest-neighbor intralayer bonds 
and the nearest-neighbor B-M bonds. This means that 
a nearest-neighbor bond model suffices to describe the 
two materials. This is not at all trivial, because in in- 
termetallic compounds, e.g., FeAL_CoAl and NiAl the 
further distant bonds are essentialliJ. 

In Table |l| the covalent bond energies for the most 
important atom pairs arc further analysed by consider- 
ing the dominant angular-resolved contributions. In all 
cases, the p-p contributions are strongest and the s-s con- 
tributions weakest. It becomes obvious from Table || 
that the d-orbitals make a non-negligible contribution to 
the bonding between the B and the M layer. Because the 
d-orbitals are not occupied in the respective free atoms 
we would denote this as a "metallic contribution" accord- 
ing to our working hypothesis of section ||. This metallic 
contribution is stronger for AIB2 than for MgB2. Never- 
theless, the "covalent contribution" is dominant also for 
the bonding between the layers, mediated by the metal 
atom. 

The energy-resolved covalent bond energies for the p-p 
and s-s contributions of the B-B intra, B-M and M-M in- 
tra atom pairs are shown in Fig. |l|, for MgB2 and AIB2. 
The benefit of the energy-resolved representation is that 
we can clearly discriminate between bonding and anti- 
bonding states. The s-s bonds are weakened because both 
bonding and antibonding states are occupied. Without 
the energy-resolved analysis we could erroneously assume 
that the s-s covalent bond energy is low because the corre- 
sponding matrix elements are small. For the other orbital 
pairs of Table [Tl| dominantly bonding states are occupied, 
but there are slight differences between MgB2 and AIB2. 
For instance, the stronger B-B intra bond energy of MgB2 
(Table |) results mainly from a stronger p — p contribu- 
tion (Table ||), and the reason for this is that for the 
trivalent Al part of the antibonding p-p states are occu- 
pied (Fig. |l]). Furthermore, the stronger B-M bond for 
M=A1 (Table again results mainly from a stronger p— p 
contribution, and the reason for this is that the bonding 
p-p states are more filled for Al than for Mg. 
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V. CONCLUSION 

We have represented a new variant of an energy- 
partitioning scheme for the density-functional total en- 
ergy which allows to define a covalent bond energy which 
is invariant against a constant potential shift. This prop- 
erty is a precondition for the comparison of the bonding 
properties in different systems within the framework of 
band-structure calculations. 

The atom-, orbital- and energy resolved contributions 
to the covalent bond energy have been calculated for 
MgB2 and AIB2. One central result of the calculations 
is that a nearest-neighbor bond model is sufficient to 
describe the bonding properties of these materials. A 
working hypothesis has been introduced to discriminate 
between a covalent and a metallic bonding character. It 
has been shown that both in MgB2 and AIB2 there is a 
metallic contribution to the bonding between the layers 
which is stronger for Al than for Mg. The benefits of the 
energy-resolved representation are demonstrated by dis- 
cussing the differences between the trivalent Al and the 
divalent Mg. 
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B-B 

intra 


B-B 
inter 


M-M 
intra 


M-M 
inter 


B-M 


total 


MgB 2 
AIB2 


-193.5 
-180.3 


1.3 
3.7 


-37.6 
-41.7 


-2.7 
10.2 


-73.7 
-86.9 


-3031 
-3179 



TABLE I. 
atom pairs. 



Covalent bond energies (in mRy) for various 





B-B intra 




B-M 




M-M intra 




p-p p-s s-s 


p-p 


p-d p-s 


s-s 


P-P 


MgB 2 


-84.9 -42.8 -14.4 


-24.9 


-13.2 -12.8 


-3.6 


-10.0 


A1B 2 


-72.9 -39.5 -14.2 


-30.8 


-17.6 -14.1 


-1.7 


-12.3 



TABLE II. The dominant orbital-resolved bond energies 
(in mRy) for various atom pairs. 
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Energy [eV] 

FIG. 1. Energy resolved covalent bond energies for the p-p (full line) and s-s (dashed line) contributions of various atom 
pairs in MgB2 (upper panels) and AIB2 (lower panels). 



G 



